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Abstract 

We present new measurements of heavy cosmic-ray nuclei at high energies per- 
formed during the first flight of the balloon-borne cosmic-ray experiment CREAM 
(Cosmic- Ray Energetics And Mass). This instrument uses multiple charge detectors 
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and a transition radiation detector to provide the first high accuracy measurements 
of the relative abundances of elements from boron to oxygen up to energies around 
1 TeV/n. The data agree with previous measurements at lower energies and show a 
relatively steep decline (~ E~°' & to E~ ' 5 ) at high energies. They further show the 
source abundance of nitrogen relative to oxygen is ~ 10% in the TeV/n region. 

Key words: cosmic rays, cosmic-ray propagation, 
PACS: 96.50.S 



1 Introduction 

The first flight of the Cosmic-Ray Energetics and Mass (CREAM) high- 
altitude balloon experiment took place in Antarctica during December 2004 
and January 2005. A central goal of this flight was to investigate the source 
properties of cosmic rays through the careful measurement of those cosmic 
nuclei which are predominantly of secondary origin. 

The escape of particles from the Galaxy during propagation is known to be 
energy dependent at high energy. As a result, the spectrum of cosmic rays mea- 
sured at Earth is different from the spectrum of these particles at their source. 
Understanding the magnitude of this difference is crucial for solving the puz- 
zle of cosmic-ray origins. Secondary nuclei are particularly useful in addressing 
this goal because they are produced largely by the spallation interactions of 
primary particles as they propagate from their source regions through the in- 
terstellar medium to Earth. The pattern of relative abundances of secondary 
and primary nuclei therefore encodes information about the propagation his- 
tory of the primary cosmic rays, allowing us to identify the contribution from 
propagation to the spectrum observed at earth. A recent review has identified 
the measurement of secondary nuclei abundances as a key goal for achieving 
future progress in cosmic-ray science p]. 

It is known from earlier measurements (see, e.g., [2J3i4j5|) that the diffusion 
escape time for particles from our Galaxy decreases with increasing particle 
energy, or magnetic rigidity. This decrease can be described as a power law 
above some threshold rigidity near 2-10 GV. Below this rigidity, the escape 
time is roughly constant and the interstellar pathlength is approximately pro- 
portional to the particle velocity /3(=v/c)[6]. Above 10 GV, the escape time 
can be characterized as a pathlength (in g/cm 2 ) by assuming a constant den- 
sity for the interstellar medium. A typical form for the rigidity dependence of 
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this quantity is A = Xq(R/ Rq)~ s ■ Here A is the mean escape pathlength, R 
is the nucleus magnetic rigidity, and Ao(~10g/cm 2 ) is the pathlength at the 
threshold rigidity R (typically 10 GV). 



The variable 5 is the key parameter in the description of cosmic-ray diffusion 
from the Galaxy at high energies. In the simplest models of cosmic-ray trans- 
port, primary particles from a source with spectral index a escape more easily 
at high energy to produce an steepening in the slope of the energy spectrum 
observed at Earth to an index of a + 5. 



The propagation of cosmic rays in our Galaxy is an essentially statistical 
process at energies where the radius of gyration in typical bulk fields is sig- 
nificantly less than the size of randomness in the magnetic field structures. 
In practice, this means that at the energies where the measurements of sec- 
ondary nuclei have thus far been possible (<100 GeV/n) the propagation is 
indeed well-approximated by models where cosmic rays simply diffuse out of 
the Galaxy. As a result, the pathlength distribution for each individual nuclear 
species is close to an exponential in form, with a mean value at high energy 
given by the expression above. 



The standard procedure when measuring secondary nuclei is to compare their 
abundances to those of their most likely (primary) parent nuclei, and examine 
how these ratios vary with energy. This can determine the value of the escape 
energy dependence, 8. For example, the ratio of boron (essentially all of which 
is secondary) to its most likely parent nucleus, carbon, is expected to decline 
with energy with a slope which almost directly yields a value for 5 above ~ 30 
GeV/n. nitrogen, on the other hand, is only partly secondary; it is known 
to have a residual abundance at the source. The ratio of nitrogen to oxygen 
is expected to decline less rapidly with energy. Finally, the ratio of the two 
primary nuclei carbon and oxygen is expected to remain relatively constant 
with energy. 



In this work, we describe measurements of the ratios of boron to carbon, 
nitrogen to oxygen, and carbon to oxygen over an energy range of lGeV/n- 
lTeV/n. Heavier nuclei (Z>8) will be the subject of future publications. One 
crucial aspect of these measurements is the correction for secondary particle 
production in the local atmospheric overburden above the balloon instrument. 
The systematic error associated with this correction becomes significant at 
high energies and ultimately imposes a limit for the maximum energy at which 
such measurements might be pursued on a balloon flight. 
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Fig. 1. A Schematic view of the CREAM instrument for the first flight. 
2 CREAM Instrument and Flight 

The CREAM instrument was first launched from Williams Field near Mc- 
Murdo Station in Antarctica on December 16th, 2004. The payload circled the 
South Pole three times in a flight which lasted nearly 42 days [7]. The mean al- 
titude during the flight was 38.6km, corresponding to an average atmospheric 
overburden of 3.9 g/cm 2 . The data for the present study were provided by a 
trigger scheme designed to tag heavy nucleus events (charge Z>3). We term 
these triggers "Hi-Z" events, around 15 million of these events are used in the 
present analysis. 

The instrument configuration used for this flight is shown schematically in 
Figure [TJ A much more detailed description of the instrument, along with a 
discussion of calibration processes, detector performance, and general archi- 
tecture is given in a recent instrument paper [8]. Here we briefly review only 
those detector elements used for the Hi-Z analysis. 

The detector system consists of a stack of horizontal components resting on a 
solid base plate. Importantly, there is no pressure gondola shell for CREAM. 
Not only does this significantly decrease the mass of the payload, it reduces 
the material overburden of the instrument, limiting the probability of local 
interactions of cosmic-ray nuclei. Such interactions can cause backgrounds of 
secondary particles via charge- changing nuclear processes. 

At the base of the detector stack is a hadronic calorimeter (CAL) which has 
been described in detail elsewhere j9]. The calorimeter is designed to measure 
particles at the highest energies and is not used in the present analysis. Above 
the CAL is a silicon charge detector (SCD) consisting of 2912 individual pixels, 
each with an area of 1.55 x 1.37 cm 2 [10]. The pixels overlap slightly to provide 
full coverage over a total area of 78 x 79 cm 2 . The dynamic range of the readout 
system can accommodate nuclei from Z = 1 to Z — 28. Over this range the 
detectors have a remarkably linear response. 



Above the SCD lies a precision gas transition-radiation detector (TRD)[TT]. 
split into two segments. Each segment consists of blocks of plastic foam (120 x 
120 cm 2 in area) in which are embedded 2 cm diameter thin-walled propor- 
tional tubes. There are 256 such tubes in a segment, each filled with a 95% 
xenon 5% methane gas mixture. The CREAM coordinate system is defined 
with the z direction vertical in Figure [H with perpendicular x and y direc- 
tions corresponding to the major axes of the detector planes. The TRD tubes 
in each segment are divided into 4 sets of 64, each of which is alternatively 
aligned with either the x or y axis. This provides a view of the particle track 
in both the xz plane and the yz plane. 



Each proportional tube in the array is instrumented with a dual-gain amplifier 
system which provides a dynamic range of ~4000. This allows for the mea- 
surement of nuclei from Z = 3 to Z = 26. Between the two TRD segments 
is an acrylic Cherenkov detector (CD), which is instrumented by eight pho- 
tomultiplier tubes (PMTs) viewing wave-shifter bars along the edges of the 
plastic. The central radiator has a refractive index near 1.5 and is doped with 
a wave-shifting dye to provide near-isotropic emission within the volume. 



At the top of the instrument is a plastic scintillation assembly, the Timing 
Charge Detector (TCD). The detector consists of eight 120 x 30 x 0.5cm 3 
paddles, viewed at each end by fast photomultipliers. They are arranged in two 
orthogonal layers of four paddles each, providing a coverage area of 120 x 120 
cm 2 . Both the timing and amplitude of the scintillator signals from each event 
are digitized by these detectors. In this work, however, we only use the event 
amplitudes. For these, the TCD has a dynamic range for nuclei from Z — 1 
to Z = 28. The total grammage of the detector stack, from TCD to SCD is 
about 7.0 g/cm 2 . 



The nuclei presented in this paper are almost all in the energy range below 
1 TeV/n. The energy measurement device used in this regime is the TRD. 
Although the Calorimeter trigger threshold accepts events near the upper end 
of this energy regime, it is not included in the present analysis since it restricts 
the available collection aperture. While the amount of transition radiation 
produced at this energy is not large (the nominal threshold for production is 
around Lorentz factors 7 > 10 3 ), the proportional tube gas, xenon, has one 
of the largest known relativistic rises in ionization energy loss[12j. This allows 
for an estimate of particle 7 in the range 4 < 7 < 5 x 10 2 . Importantly, the 
response of this detector has been calibrated during beam tests at CERN. 
This is discussed extensively elsewhere [T3"] . 
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Fig. 2. Scatter plots of detector signals for a small fraction of the flight data. Left: 
The square root of the TCD signal versus the square root of the CD signal. Also 
shown are representative cuts applied to select relativistic nuclei (between dashed 
lines). The vertical scale here has been adjusted to approximate the nuclear charge 
and all signals are corrected for gain and pathlength. Right: The square root of the 
corrected TCD signal versus the square root of the corrected SCD for the selected 
relativistic nuclei from the left panel. Again, representative cuts are shown which 
select nuclei that do not undergo charge-changing interactions while traversing the 
instrument. 

3 Data Overview 



In this section we describe in some detail the data collected during the flight, as 
well as some of the corrections applied to those data. A Hi-Z event is required 
to generate a signal above a fixed threshold in both the CD and TCD within 
a coincidence window of ~ 100ns. These thresholds were adjusted before the 
flight to ensure high efficiency for relativistic boron [Z = 5) nuclei, while 
rejecting the tail of the much more abundant Helium nuclei (Z = 2). 



3.1 Detector Corrections 



One of the benefits of a tracking instrument is that it provides a simple and 
direct mechanism for identifying and correcting non-uniformities in detector 
response. After constructing particle trajectories using the tracking algorithm 
described below, we apply mapping corrections for every event collected. The 
maximum level of these corrections is ~50% for the TCD and ~30% for the 
CD. No positional mapping is required for the silicon detector, but the in- 
dividual pixels do receive a gain-matching correction. Again using tracking 
information, the signals from all the detector systems are also corrected for 
pathlength. 
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In the TCD, an additional correction is made for the (Birk's Law) saturation 
in the scintillator and non-linearity of the PMT response. The TRD's propor- 
tional tube signals are also corrected for gain variations caused by fluctuations 
in pressure and temperature. These signals are gain-matched to a precision of 
~3% over the entire flight. Additionally, there were several calibration proce- 
dures implemented periodically throughout the flight, as discussed in [8]. 



3.2 Charge Identification 

To outline the analysis procedures, we first discuss in general terms the identi- 
fication of charges using the CREAM detectors. This initial estimate of charge 
is used as an input to one of the tracking algorithms. Later, in Section HI we 
describe the full relativistic energy analysis which requires more detailed inter- 
pretation of the events. The left-hand panel of Figure [2] shows a small sample 
of events measured in the TCD and CD after applying the above corrections. 
The individual charges are clearly visible. The vertical TCD scale is arbitrary, 
but has been adjusted to approximate the derived charge of the incoming 
relativistic nuclei. Also apparent on this plot is the threshold setting for the 
CD (abscissa), which can be seen to admit relativistic boron nuclei, and with 
lesser efficiency, Beryllium. 

In a polar flight, the geomagnetic cutoff rigidity is very low, admitting a large 
flux of low-energy particles. Since we are primarily interested in high-energy 
particles for this study, we use the Cherenkov Detector to identify and reject 
the low energy ones. The CD has a refractive index of 1.5, which provides a 
threshold velocity of (3 ~ 0.7. Particles with low CD signal are therefore easily 
identified as having an energy of a few hundred MeV/n, well below minimum 
ionization in the proportional tubes of the TRD. An example of the cuts used 
to reject such nuclei are shown in the left panel of Figure [2] (dashed lines). 

Also shown in Figure [2J (right panel) is a scatter plot between the TCD and 
the SCD signals (after the above cuts are applied). In general, the charge 
measurements by these devices are in good agreement. However, there are a 
number of events in which the charge measured by the SCD (abscissa) is lower 
than that measured by the TCD. These are caused by nuclear interactions in 
the instrument. Such events represent a small fraction of the data [14] and can 
be easily removed from the data by a simple consistency check between the 
TCD and SCD signal. A representative cut is shown by the dashed lines on 
the figure. 

After the relativistic particle cuts are applied to the CD data and the con- 
sistency cut has been applied between the TCD and SCD signals, the overall 
charge resolution can be determined. Figure [3] shows the mean of the SCD 
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Fig. 3. Charge histogram of nuclei with Hi-Z triggers. The charge estimate shown 
is the mean of the square root of the corrected TCD signal and the square root of 
the corrected SCD signal. The data were selected using the cuts from Figure [2j The 
left-hand panel has a linear scale and the right-hand shows the same data with a 
logarithmic scale. The distribution for carbon has a width of <r^=0.16 charge units. 

and TCD signals ((y/TCD + \fSCD)/2) for relativistic Hi-Z, non-interacting 
particles collected during the flight. The charge distribution extends down to 
Lithium, which produces Hi-Z triggers at low efficiency because of the fixed 
CD threshold discussed above. Across the whole charge spectrum, all of the 
nuclei are visibly separated from their neighbors, including the rarer nuclei, 
Fluorine (Z=9) and Sodium (Z=ll). Overall, the achieved charge resolution 
is excellent, with an RMS (la) resolution for carbon nuclei of AZ = 0.16 
charge units. Performance at this level is crucial for the successful study of 
rare secondary particles, as discussed below. 



4 Analysis Methods 

In the following section we will discuss the analysis procedures and corrections 
used to reduce the data from the CREAM flight. This includes a discussion 
of the final selection cuts (which are somewhat more restrictive than the cuts 
discussed above), the energy-determination mechanism, and energy deconvo- 
lution procedures. We begin with a synopsis of the track-fitting algorithms. 

4 . 1 Event Track Fitting 

The heart of the analysis is the track fit. The results of this fit are used in the 
mapping corrections described above to correct detector signals. They are also 
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used to determine the actual track length of the particle in each of the TRD 
tubes it traverses. The accuracy of this information is crucial for properly 
determining the dE/dx signal for each event. Because this number needs to 
be determined as well as possible, the fitted tracks must have a precision 
significantly less than a tube diameter. This is achieved by a multistage track- 
fitting procedure. 

The first-stage track fitter employs a simple and very fast least-squares mini- 
mization which ignores the physical extent of the proportional tubes. Outlying 
signals are rejected from the fit by iteratively contracting an envelope around 
the estimated track position. Within 3 passes, the precision of this initial fit 
is ~5mm FWHM in both the xy and yz planes. The fit is then improved 
using a more sophisticated second-level algorithm which accounts for the full 
three-dimensional structure of the TRD. This algorithm uses the Minuit soft- 
ware package to perform a maximum-likelihood fit incorporating the different 
energy-loss distributions produced by tracks at different impact parameters. 
Because these distributions also vary with particle charge, an initial determi- 
nation of particle charge Z (discussed above) is used as input to the algorithm. 
When initialized with the parameters provided by the level-one track, this fit- 
ter achieves an accuracy at the level of ~1 mm RMS. 

Using the tracking information, the TRD's dE/dx signal is directly derived 
from the individual tube signals. A quality cut is applied which requires in- 
dividual track segments to be at least 1.5 cm long, helping to minimize the 
impact of fluctuations for thin gas layers. A simple test of the accuracy of this 
process can be obtained by comparing the signal derived from the xz plane to 
that of the yz plane. These fits are essentially independent except for the com- 
mon track energy-loss parameter, and therefore a comparison gives a rough 
estimate of the overall error in the determination of the dE/dx value. The left 
panel of Figure H] shows such a test for flight oxygen nuclei. The relative RMS 
difference between the xz and yz determinations is ~ 14%, implying that the 
resolution in the combined TRD is a ~7%. The actual resolution, discussed 
below, is a little broader than this because of the correlations introduced by 
the track energy-loss parameter. 



4-2 Energy Determination 

The energy of incoming particles can be determined in two ways. At lower 
energies, the signal from the Cherenkov Detector can be used to estimate 
velocity, and hence energy. The TRD proportional tubes are used at energies 
above where the CD response has saturated. The signal in these tubes is 
produced by a combination of two energy-loss mechanisms, both of which scale 
with the energy of the nucleus. The primary component comes from Coulomb 
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Fig. 4. Left: The fractional difference between the TRD signal determined from the 
xz track and the yz track for oxygen nuclei. The result has an RMS value of ~14%. 
Right: The TRD sum signal distribution for flight oxygen nuclei, selected to be at 
minimum ionization using the CD. The RMS resolution is a ~10%. 



10 



collisions of the nuclei with electrons in the gas. An additional contribution 
comes from any x-ray transition radiation absorbed in the xenon gas mixture. 
In the energy region of the present measurements (< 1 TeV/n), the tube 
signals are essentially all due to ionization energy loss. 

The right panel of Figure H] shows the distribution of TRD signals for oxygen 
nuclei which are selected to be close to minimum ionizing (using the Cherenkov 
detector) . These events are also selected to have more than six TRD tubes with 
a usable signal. It is the width of this distribution (which remains essentially 
constant throughout the relativistic-rise region) that ultimately determines 
the energy resolution of the TRD. The RMS (la) width of the signal, S, 
distribution is AS ~ 8 in these units. As demonstrated in Figure [51 which 
shows the response function of the TRD, this corresponds to a resolution in 
logio(7) ~ 0.3, where 7 is the particle Lorentz Factor. This figure also shows 
the good agreement of the GEANT4[15] detector simulation (small points) 
with the results of a test beam calibration of the same configuration (large 
points). This calibration is discussed in more detail elsewhere [T3f8] , and has 
been shown to be consistent with the energies derived from the calorimeter 
signals in the ~100 GeV/n region for carbon and oxygen[16j. 

In practice, a simple least-squares energy fitter is used to determine the particle 
energy which best corresponds to the signals measured in the TCD, CD and 
TRD. For each event, a goodness-of-fit, x%> * s generated, based on the known 
response functions of these detectors. 



4-3 Selection Cuts 

The data for this study were selected from time periods when the experiment 
was relatively stable and when there was no recent gas venting, filling or other 
operations, and the atmospheric overburden was <5g/cm 2 . Periods were also 
avoided when the high voltage was being adjusted on any experimental system. 
These selections resulted in an effective data collection time of ~21 days. 

The events used in the final analysis must pass various identification and 
quality cuts. The nuclear charge is selected by placing cuts on the TCD and 
SCD signal distributions, similar to those shown in Figure [2j An additional 
tracking cut helps here by requiring that the calculated particle trajectory 
pass through the SCD within 2 cm of the pixel with the largest signal. This 
cut keeps only "clean" events and eliminates possible albedo effects from the 
calorimeter just below the SCD. Approximately 70% of Hi-Z triggers with 
TRD data survive this cut. To provide adequate signal resolution, only events 
with more than six analog samples in the tube array are selected for analysis. 
In addition, the total track pathlength in the tubes is required to be at least 
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Fig. 5. The response of the TRD vs particle log(Lorentz factor). The small sym- 
bols are from a GEANT4 simulation of the detector configuration, and the large 
symbols are data from a direct accelerator measurement, as discussed in [13J. Also 
shown schematically is the relationship between the TRD signal resolution and the 
log(Lorentz factor). 
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12 cm. These two cuts admit ~25% of the events with a good level-two track 
fit. A consistency check between the xz-pl&ne TRD signal and the yz-view 
TRD signal (amounting to a simple cut on the distribution shown in Figure H]) 
further improves the overall resolution, while accepting ~65% of the remaining 
events. A final cut is made on the likelihood of the track fit, rejecting poor 
geometry events, while keeping over 90% of those events which make it to this 
stage. The cumulative efficiency of all the cuts applied to the events which 
trigger the instrument in this analysis is ~10%. 

As mentioned above, the energy fitter produces a simple Xe value which can 
be used to reject pathological events. A loose selection cut is made that keeps 
only those events with reasonable values of x\ < 10- This works rather well 
except when knock-on electrons, or other processes, produce non-statistical 
upwards fluctuations in the CD signals. This effect can cause a nucleus at low 
energy to mimic a high-energy event in the TRD. This is easily understood 
by referring to the TRD response curve shown in Figure [51 There is clear 
ambiguity in energy for a TRD signal of, for example, 115 scale units. The 
CD is needed to determine on which side of minimum ionization the event lies. 
A cut is therefore placed on the CD to ensure the events are on the high-energy 
side. 

The flight data themselves can be used to estimate the effectiveness of this 
Cherenkov cut. In Figure E] the CD signal is plotted for oxygen nuclei which 
are selected as discussed above except for the Cherenkov cut illustrated in 
Figure [2j These also are required to have a TRD signal between 30 and 35% 
above minimum ionization. According to Figure these could either have 
an energy ~400 MeV/n or ~100 GeV/n from the response curve. It can be 
readily seen that there are two easily separable populations - corresponding to 
the low-energy and high-energy events. A simple cut on the CD signal (above 
1.05 in this case) ensures only the high-energy particle survive. 



4-4 Charge Overlap Rejection 

Reliably discriminating between different charges is extremely important for 
any experiment measuring secondary nuclei. For example, because boron at 
high energy has only a few percent of the abundance of carbon, the peaks of 
the reconstructed populations of these nuclei must be at least three standard 
deviations apart to prevent a significant background from carbon from leaking 
into the boron sample. 

With a TRD instrument, this requirement becomes even more stringent be- 
cause the TRD signal response is not independent of Z (it scales as Z 2 ). As a 
result, if a low-energy carbon nucleus gets misidentified by the charge detec- 
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Cerenkov Signal Size TCD Signal Size 

Fig. 6. Left: The distribution of Cerenkov signals from oxygen events selected to 
be between 115 and 125 on the response curve of Figure See the text for details. 
The right-hand population of events comprises particles near 100 GeV/n and the 
left-hand peak is a tail of low-energy events near 400 MeV/n which have passed 
identical cuts. A cut at 1.05 provides effective discrimination between the two pop- 
ulations. Right: The distribution of charge signals measured by the TCD for events 
which have an SCD charge of Z = 5.0 ± 0.2. The effects of spallation are obvious as 
several of the nuclei enter the instrument with a charge larger than that measured 
in the SCD. 
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tors as a boron, it will not be counted as a low-energy boron, but rather as a 
high-energy boron because of the incorrect Z 2 normalization factor. This exac- 
erbates the "charge overlap background" problem because the steeply falling 
cosmic-ray spectrum means there are thousands more low-energy carbon nu- 
clei than high-energy boron nuclei. 

Due to this additional effect, a separation of at least four standard deviations 
between the boron and carbon populations is required to make a background- 
free determination of the boron abundance. A similar condition applies to the 
measurements of nitrogen and oxygen. 

The CREAM instrument fortunately has two independent charge-determining 
detectors, the TCD and SCD, each of which has a resolution significantly 
less than 0.2 charge units (over the range of charges considered here). To 
combat the charge overlap problem, we place circular cuts of radius AZ=0.2 
on all the charge peaks of the right panel of Figure [21 This selection produces 
effective cuts at AZ=±0.14 on the summed charge distribution and completely 
removes interacting events. The resulting separation between the cut charge 
populations exceeds 4.5 standard deviations, thereby providing a clean sample 
of the low-abundance boron and nitrogen nuclei at high energies. 

4-5 Energy Deconvolution 

All the events passing the above cuts are collected into energy bins of width 
logio(E)=0.4. This corresponds to a width of ~ 1.3<r in terms of the energy 
resolution, discussed above. Approximately 45,000 relativistic carbon nuclei 
and 40,000 oxygen nuclei make it to this stage. 

The true energy spectrum is derived from these binned values by correcting 
for energy spill-over between the bins due to the finite instrumental resolution. 
To account for these effects in a system where the resolution is a function of 
energy, a simulation is generally required. For this measurement, the primary 
energy determination device is the TRD, so we use an accurate simulation 
of the TRD and its fluctuations. This simulation uses a GEANT4 model of 
the upper CREAM payload including the TCD, TRD and SCD. The TRD 
component of this simulation has been discussed previously in [T3] . 

The events from this simulation are passed through the same analysis code 
that is used to evaluate and fit the flight data. Various cross-checks are made to 
ensure the simulation provides an accurate representation of the data produced 
in the flight, including both mean signal levels and the fluctuations of those 
signals. Using the simulated results, a transfer matrix is obtained which can 
be used to deconvolve the measured energy bins to produce estimated true 
energy bin populations. This method, which has been discussed in previous 
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Fig. 7. Measurements of the ratios of nuclei as a function of energy. Left: Filled circles 
show the ratio of boron to carbon vs. energy after corrections. The horizontal errors 
are an estimate of the systematic error in the overall energy scale. The thin vertical 
lines correspond to the statistical error of the ratio and the grey bars show the 
systematic uncertainty in the ratio. See text for details. The lines represent model 
calculations for various values of the magnetic-rigidity dependence parameter, 5, 
in escape from the Galaxy - as discussed in the text. These are; solid line 6=0.6, 
long-dashed line J=0.333, short-dashed line 5=0.7. The stars are data from the space 
experiment, HEAO-3-C2[5j. Right: Filled circles show the ratio of nitrogen to oxygen 
vs. energy after corrections. The error bars and data points are as in the left panel. 
The lines represent model calculations of this ratio with the escape parameter 5=0.6 
(solid line in the top left-hand panel). The different curves correspond to different 
assumptions on the amount of nitrogen in the source material. These are; solid line 
source N/O = 10%, long-dashed line source N/O = 5%, short-dashed line N/O = 
15%. 



studies (e.g., |4j), is an iterative process since, for accurate results, the spectral 
index used in the simulation must match the (unknown) index of the data. 
However, as noted by several previous authors [HUE], the convergence of the 
iteration is rapid and small residual discrepancies have little impact on the 
final result. 
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Fig. 8. Left: Filled circles show the ratio of carbon to oxygen vs energy after correc- 
tions. Right: The fluxes (energy/nx area x time x solid angle) -1 (arb. units) of carbon 
nuclei (open triangles) and oxygen nuclei (open circles) from this measurement. The 
data are consistent with the dashed line, a power law oc(Energy/n) _a with a spectral 
index a=2.6. 

5 Results and Corrections 



Using the analysis, selection cuts, and deconvolution scheme outlined in the 
previous sections, the flight data are processed to determine the ratios of ele- 
mental abundances as a function of energy into the TeV/n region. As described 
above, this is achieved with a simple binning and deconvolution process which 
assigns the median energy of each bin by assuming a spectral index close 
to that of each denominator element. These results are given in Tables 1-3, 
and the associated Figures [7| and [SJ Also shown for comparison are the re- 
sults from the space measurement HEAO-3-C2[5j which used an assembly of 
Cherenkov counters of different thresholds to measure elemental composition 
up to ~ 35GeV/n. These are the most statistically significant measurements 
to date in the lower energy region. These tables show the abundance ratios 
B/C, N/O, and C/O, respectively, as measured at the instrument as well as 
those values after extrapolation to the top of the atmosphere. 

The corrections to the top of the atmosphere are made by considering sepa- 
rately the probability of charge- changing interactions in the instrument and 
the atmosphere. The right panel of Figure [6] clearly shows that events that 
have fragmented can be eliminated with charge cuts. However, the efficiency 
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of these cuts must be accounted for in the analysis procedure. 

Charge-changing interactions also occur in the atmosphere above the instru- 
ment. Obviously, these events cannot be identified, so a post-analysis correc- 
tion must be applied using nuclear interaction models. For this purpose, we 
utilize both the partial cross-section model of Tsao et al. [19] and GEANT4 
nuclear interaction simulations. These corrections become extremely impor- 
tant as the elemental abundance ratios decline. For example, if 3% of carbon 
nuclei (Z=6) change to boron (Z=5) in the atmosphere above our instrument 
this is a fairly negligible correction at 5 GeV/n, where B/C~30%. At around 
~ ITeV/n, however, it becomes significant, since in this region B/C~5%. 

Since it appears that the B/C ratio continues to decline relatively steeply 
with energy, the accurate determination of this ratio by balloon payloads will 
be intrinsically limited to the region below a few hundred GeV/n because of 
uncertainties in these corrections. Measurements well above these energies can 
probably only be made successfully on a future space mission, where there are 
no systematic limitations from atmospheric corrections. 

There are several different uncertainties documented in the tables. For the 
energy values, the quoted errors are systematic - arising from the precision with 
which the data and simulation response curves can be aligned. We estimate 
the fractional error in this process to be ~3%. In the response region where 
most of these data were collected, this precision corresponds to an energy-scale 
uncertainty of ~15%, as displayed in the tables. 

The entries for the corrected abundance ratios feature both statistical and 
systematic errors. The statistical portion is derived in the standard way by 
propagating the statistical errors of each element in a ratio, using the original 
counts of events in the bins before deconvolution and before the atmospheric 
corrections. The systematic errors in each ratio stem from two effects. The 
first stems from residual uncertainty in the charge- dependent efficiency of the 
selection cuts used to construct the ratios. There are several cuts for which 
the charge dependence may amount to a few percent, so we use a conservative 
value of 10% for the total charge- dependent efficiency error. A second system- 
atic effect derives from uncertainties in the atmospheric secondary corrections. 
The partial cross-section values in the literature for carbon+Air — > boron in 
the relativistic region can differ by up to ~30%. The typical correction in the 
ratio of B/C from atmospheric secondaries is of the order ~ 3%, and we as- 
sign a fixed 1% systematic error for the uncertainty in correcting the ratios to 
the top of the atmosphere. The systematic errors quoted in the tables reflect 
whichever contribution dominates in the appropriate energy region. 

To provide a check on these measurements we also present data on the C/O 
ratio in Figure [SJ As both C and O are primary nuclei - originating mostly in 
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the cosmic-ray source - the expectation is that their abundance ratio should 
be close to constant with energy. This expectation is in agreement with our 
measurements. 

A further check can be seen in the right panel in Figure El Here the measured 
intensities of carbon and oxygen nuclei are shown as a function of energy (these 
are differential elemental intensities given in arbitrary units since the overall 
live-time of the flight has not yet been accurately determined). Apart from 
demonstrating that the CREAM measurements reported here cover nearly 9 
decades in differential intensity these data show that the measured energy 
spectra for primary nuclei are consistent with a differential energy spectral 
index of a=2.6 at energies above lOGeV/n. This value is near the average of 
previously-reported measurements and provides confidence that our analysis 
techniques and energy scales are accurate. 



6 Discussion 

The left panel of Figure [7] shows that the measured ratio of B/C declines fairly 
steeply with energy into the TeV/n region. Overlaid on this figure are three 
lines which represent the results of a simple leaky-box model of cosmic-ray 
propagation. For each of the lines, a different value of 5 is used. The upper- 
most (long-dashed) line uses 5=0.33, the solid line uses 5=0.6, and the short- 
dashed line uses 5=0.7. The small data points (stars) on the figure represent 
the ratio measurements of the HEAO-C2-3 experiment [5] . The present data 
are consistent with these measurements and also seem to be consistent with 
the 5=0.6 curve. 

Measurements of N/O allow one to investigate the apparent source abun- 
dance of nitrogen in cosmic rays, a topic which has been discussed for many 
years. The cosmic nitrogen abundance is somewhat unusual in that it con- 
tains significant contributions from both primary and secondary sources. At a 
high-enough energy the primary component should come to dominate and the 
ratio N/O should flatten out. Our N/O measurements, again together with 
those of HEAO-C2-3 are shown in the right panel of Figured As before, three 
simple leaky-box propagation models are overlaid on the data, in this case 
with the parameter 5 fixed at 0.6. The different lines on this plot correspond 
to different assumptions on the source abundance of nitrogen compared to 
oxygen. The solid line corresponds to a source N/O = 10%, the long-dashed 
line corresponds to a source N/0=5% and the short-dashed line is for source 
N/O = 15%. Although lower-energy measurements on nitrogen isotopes have 
seemed to favor a source abundance N/O ~5% (e.g., [20]), these new data, 
which agree well with HEAOJS], seem to prefer a source value ~10%. 
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7 Conclusions 



We have presented new measurements of the relative abundances of cosmic-ray 
secondary nuclei at high energies, where the energy dependence of cosmic-ray 
propagation is relatively unexplored experimentally. These results have been 
produced with the first flight of the CREAM high-altitude balloon instru- 
ment in Antarctica during austral summer 2004/2005. These are challenging 
measurements, and the results span over three decades in energy and over 
nine decades in particle intensity. The energy scale for these results is pro- 
vided by a transition-radiation detector, operating mostly in the regime of the 
relativistic rise. A critical part of the success of this effort is the excellent res- 
olution of the CREAM charge detectors. Both the TCD (plastic scintillator) 
and SCD (silicon detector) have demonstrated charge resolutions of AZ<0.2 
charge units. This provides excellent charge isolation between primary nuclei 
and the much less abundant secondary nuclei at all energies. Indeed, without 
this level of charge separation, the background of misidentified nuclei would 
make this measurement difficult if not impossible above ~100 GeV/n. 

The results presented here show the interstellar propagation pathlength de- 
creases fairly rapidly with energy, with an energy dependence in the range 
5 ~0.5-0.6. As a result, the propagation pathlength of cosmic nuclei is smaller 
by an order of magnitude for particles in the TeV/n region compared to those 
at lower energies below 10 GeV/n. This high-energy pathlength (~1 g/cm 2 ), 
while small, is still large compared to the typical grammage of the Galactic 
disk (< 0.002 g/cm 2 ), and so it does not significantly constrain residual path- 
length models proposed for higher energies [2 lj or in the source region [22J. The 
observations of N/O abundance ratios are consistent with this pathlength de- 
pendence and they suggest an N/O source abundance close to 10% - larger 
than some previous estimates based on lower-energy isotope measurements 
(see e.g. [23]). 

For this measurement, the dominant source of systematic error at high energy 
are uncertainties in the cross sections for producing secondaries by charge- 
changing interactions in the atmosphere above the instrument. Because there 
is a significant decrease in the interstellar pathlength with energy, the amount 
of boron, for example, at high energy is small - making the impact of uncer- 
tainties in the atmospheric boron contribution significant above ~100 GeV/n. 
At these energies the contribution from charge-changing interactions in the 
atmosphere is similar in size to the total production of boron during prop- 
agation through the Galaxy. As a consequence it seems likely that accurate 
measurements of B/C on high-altitude balloon experiments will be ultimately 
limited by systematic errors in the TeV/n region until knowledge of these 
atmospheric corrections can be improved. 
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Table 1 

Table of measurements of Boron to Carbon Ratio 



Kinetic Energy 


Ratio h>/C 


T~ ) i • , -T> / /~1 

Ratio B/C 


Energy Range 


(GeV/n) 


measured 


corrected 


(GeV/n) 


1.4±0.2(sys.) 


0.353 


0.320±0.007(siai.) ± 0.03(sys.) 


1-4 


5.7±0.9(sys.) 


0.256 


0.225±0.004(stat.) ± 0.02(sys.) 


4-16 


23±3(sys.) 


0.185 


0.155±0.005(stat.) ± 0.014(sys.) 


16-63 


91±14(sys.) 


0.131 


0.101±0.011(stot.) ±0.01(sys.) 


63-251 


363±54(sys.) 


0.099 


0.071±0.025(stat.) ± 0.01(sys.) 


251-1000 


1450±217(sys.) 


0.061 


0.033±0.025(stat.) ± 0.01(sys.) 


1000-4000 


Table 2 








Table of measurements of Nitrogen to Oxygen Ratio 




Kinetic Energy 


Ratio N/O 


Ratio N/O 


Energy Range 


(GeV/n) 


measured 


corrected 


(GeV/n) 


1.4±0.2(sys.) 


0.341 


0.299±0.006(stat.) ± 0.03(sys.) 


1-4 


5.7±0.9(sys.) 


0.286 


0.246±0.004(stat.) ± 0.025(sys.) 


4-16 


23±3(sys.) 


0.248 


0.210±0.009(stat.) ± 0.02(sys.) 


16-63 


91±14(sys.) 


0.211 


0.174±0.026(stat.) ± 0.02(sys.) 


63-251 


363±54(sys.) 


0.160 


0.124±0.072(stat.) ±0.01(sys.) 


251-1000 


1450±217(sys.) 


0.083 


0.050±0.034(stat.) ± 0.01(sys.) 


1000-4000 



Table 3 

Table of measurements of Carbon to Oxygen Ratio 



Kinetic Energy 


Ratio C/O 


Ratio C/O 


Energy Range 


(GeV/n) 


measured 


corrected 


(GeV/n) 


1.4±0.2(sya.) 


1.19 


1.10±0.01(stai.) ± 0.1(sys.) 


1-4 


5.7±0.9(sys.) 


1.19 


l.ll±0.01(stai.) ± 0.1(sys.) 


4-16 


23±3(sys.) 


1.25 


1.16±0.03(stat.) ± 0.1(sys.) 


16-63 


9l±U(sys.) 


1.34 


1.25±0.10(stai.) ±0.1(sys.) 


63-251 


363±54(sys.) 


1.35 


1.25±0.32(stai.) ± 0.1(sys.) 


251-1000 


1450±217(sys.) 


0.71 


0.66±0.41(stat.) ± 0.1(sys.) 


1000-4000 
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